Synchronization transitions on scale-free neuronal networks due to finite information transmission 

delays 
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We investigate front propagation and synchronization transitions in dependence on the information transmis- 
sion delay and coupling strength over scale-free neuronal networks with different average degrees and scaling 
exponents. As the underlying model of neuronal dynamics, we use the efficient Rulkov map with additive noise. 
We show that increasing the coupling strength enhances synchronization monotonously, whereas delay plays a 
more subtle role. In particular, we found that depending on the inherent oscillation frequency of individual neu- 
rons, regions of irregular and regular propagating excitatory fronts appear intermittently as the delay increases. 
These delay-induced synchronization transitions manifest as well-expressed minima in the measure for spatial 
synchrony, appearing at every multiple of the oscillation frequency. Larger coupling strengths or average de- 
grees can broaden the region of regular propagating fronts by a given information transmission delay and further 
improve synchronization. These results are robust against variations in system size, intensity of additive noise 
and the scaling exponent of the underlying scale-free topology. We argue that fine-tuned information trans- 
mission delays are vital for assuring optimally synchronized excitatory fronts on complex neuronal networks, 
and indeed, they should be seen as important as the coupling strength or the overall density of interneuronal 
connections. We finally discuss some biological implications of the presented results. 
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I. INTRODUCTION 

Synchronization phenomena are ubiquitous in nature and 
play an important role in biology, ecology, climatology, so- 
ciology, technology, and even fine arts 01 HI Hi- in the 
study of nonlinear dynamical systems, synchronization is re- 
currently being placed in the focus of attention, and recently 
insightful findings regarding the synchronization phenomena 
on complex networks were reported [H H Hi 1§ and 
comprehensively reviewed [3]. It is interesting to see that 
one can literary infer topological scales of complex networks 
based solely on synchronization ill ill , thus making a closed 
loop of dependence between the synchronizability and the 
structure of underlying interactions of network elements. In 
neural systems, in particular, the interplay between the net- 
work structure and the dynamics taking place on it is closely 
interrelated. The function-follow-form para digm, fo r exam- 
ple, is central to attest to this observation M12lll3lll4ll5lll6ll . 
It is also well known that the cerebral cortex features proper- 
ties that are characteristic for complex networks lfl7l flii Il9ll . 
As a results, the firing activities of individual neurons are of- 
ten related to the synchronization of the underlying network, 
and accordingly synchronized firings can be observed at vir- 
tually all processing levels, including the retina 
lateral geniculate nucleus [21], and the cortex 1122 ,1 

Synchronization on complex networks, therefore, has actu- 
ally become a focal topic in theoretical neurosciences 027112811 . 
as evidenced by several recent studies devoted to the explo- 
rations of this subject (H, HH H2, HE H3, HH ^" 
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Important works have been elaborating on gene ral aspects 
of synchronization on scale-free |T3SL l39i l40l kill as well as 
complex-gradient ll42ll networks, among many other models. 
It is now clear that synchronization is key to the efficient pro- 
cessing and transmission of information across a nervous sys- 
tem such as the brain 11231 12611 . The handling of informa- 
tion transmission over a neuronal network, however, is still 
an open avenue for research. Since information transmission 
delays are inherent to neuronal systems because of the finite 
speeds at which action potentials propagate across neuron ax- 
ons, and due to the time lapses occurring in both dendritic and 
synaptic processing [43], studies are in need of catching up 
with the most recent advances in synchronization research on 
complex networks. 

Indeed, delays have been found responsible for several in- 
teresting phenomena in coupled dynamical systems. For ex- 
ample, Ernst et al. B4411 have identified mechanisms of syn- 
chronization among pulse-coupled oscillators in the presence 
of time delay. Moreover, it has been shown that coupled os- 
cillators undergo a transition towards amplitude death faster 
if the time delays in coupling are distributed over an interval 
rather than being uniform throughout the system 14511 . The 
role of delays and connection topologies for the synchroniza- 
tion of coupled chaotic maps has been studied in [46], where 
it is reported that on scale-free and random networks suffi- 
ciently large coupling strengths can offset the delayed flow 
of information. It has also been shown that networks with 
delays can sometimes synchronize more easily than in their 
absence, and it has been argued that this may be particularly 
relevant for neuronal networks for establishing a concept of 
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collective information processing in the presence of delayed 
information transmission M47I1 . In the present paper we give 
further support to the latter assumption by considering the 
impact of delays in scale-free neuronal networks. More re- 
cently, the role of delays by the formation of the so-called 
chimera, i.e. coexistent coherent and incoherent states in a 
system of nonlocally coupled phase oscillators has also been 
examined B48I1 . and it has been shown that time-delays can 
induce a transition towards phase clustering, giving rise to 
clustered chimera states that have spatially distributed phase 
coherence separated by incoherence with adjacent coherent 
regions in antiphase. Here we support the theory of delay- 
induced dynamical transitions in terms of the synchrony of 
neuronal noise-induced excitations on a scale-free network. 
Notably, previous studies have already considered particularly 
the neuronal dynamics on large networks in conjunction with 
information transmission delay, but the focus was primarily 
on the bifurcation structure of transitions between different 
delay-induced states, including oscillatory bumps, aperiodic 
regimes, traveling, lurching and standing waves, as well as 
regimes of multistability l49ll50ll . The impact of information 
transmission delay on neuronal synchronization, on the other 
hand, has been studied in IBH l52l l53ll . where the emergence 
of zigzag fronts, clustering antiphase synchronization and in- 
phase synchronization on regular and small-world neuronal 
networks has been discussed. 

At present, we aim to extend the scope of the above- 
mentioned investigations by studying front propagation and 
synchronization transitions in dependence on the coupling 
strength and information transmission delays over scale-free 
neuronal networks with different average degrees and scal- 
ing exponents. Notably, it has been reported that, by us- 
ing functional magnetic resonance imaging, power-law distri- 
butions can be obtained upon linking correlated fMRI vox- 
els 111811 . and that the robustness against simulated lesions 
of anatomic cortical networks relies mostly on the scale-free 
structure llsill . This study thus addresses a relevant system 
setup that is still widely open for new research. More specif- 
ically, we report several non-trivial effects induced by finite 
delay lengths, and the ability of its fine-tuning towards highly 
synchronized fronts of excitations. These findings are com- 
pared to the impacts of different coupling strengths and av- 
erage degrees, and their robustness is examined at different 
levels of additive noise, variations in system size and differ- 
ent scaling exponents of the underlying scale-free topology. 
Remarkably, we found that, irrespective of the system size 
and the scaling exponent, properly adjusted information trans- 
mission delays play a pivotal role in warranting synchronized 
fronts of excitations on noisy scale-free neuronal networks, 
which can be further enhanced via larger coupling strengths 
or higher average degrees of the constitutive nodes. We argue 
that this is primarily attributed to the emergence of locking 
between the delay and the inherent oscillation frequency of 
individual neurons of the scale-free network. 

The remainder of this paper is organized as follows. In the 
next section, we describe the Rulkov map ll55ll . which will be 
employed to obtain an efficient setup for simulating neuronal 
dynamics on scale-free networks B56I1 . In Section II we also 



present the coupling scheme and the measure for synchroniza- 
tion of excitatory fronts, as well as other mathematical meth- 
ods to be used. In Section III we present the main results, 
and in the last Section we summarize our findings and discuss 
their potential implications. 

II. MATHEMATICAL MODEL AND SETUP 

For simulating the neuronal dynamics on a scale-free net- 
work effectively, the Rulkov map [55] is employed, which 
succinctly captures all the major dynamical features of the 
complex continuous-time models. The spatial-temporal evo- 
lution of the studied network, corrupted with additive Gaus- 
sian noise and experienced with information transmission de- 
lays, is described by the following iteration equations 

x {t) (n+l) = af[x®(n)] + y (i) (n) + w£ (i) {n) 

+D^e l < j [x 3 {n-T)-x l {n)] , (1) 
3 

y (i) {n + l) = y {i) {n) -px®(n) -7, i= 1,...,N 

where n is the discrete time index, a;W (n) is the membrane 
potential and j/W (n) the variation of ion concentration of the 
i-th neuron, representing the fast and the slow variable, re- 
spectively. The slow temporal evolution of y"> (n) is due to 
the small values of the positive parameters (3 and 7, which 
within this study equal p = 7 = 0.001 unless stated other- 
wise. Moreover, a is the main parameter determining the dy- 
namics of individual neurons on the scale-free network. Ac- 
cording to lf55ll . if a < 2.0 all neurons are situated in ex- 
citable steady states [x* — — l,y* = — 1 — (a/2)], whereas 
if a > 2.0 complex firing and bursting patterns of temporal 
activity emerge via a Hopf bifurcation. Here we set a = 1.95 
and initiate each neuron from steady state initial conditions, 
so that the additive spatiotemporal Gaussian noise £j(n), hav- 
ing mean < & >= and autocorrelation < £i(ri)£j(h) >— 
5ijd(n — h), acts as the only source of large-amplitude excita- 
tions. Moreover, in Eq. (1) f(x) = j^p- is a nonlinear func- 
tion warranting the essential ingredient of neuronal dynamics, 
parameter w determines the noise intensity, D is the coupling 
strength, and r is the information transmission delay. The lat- 
ter two parameters will be in the focus of attention within this 
work, whereas w, f3 and 7 will be varied only occasionally. 

As the interaction base between neurons we primarily use 
the scale-free network generated via growth and preferential 
attachment as proposed by Barabasi and Albert 115611 . consist- 
ing of N = 200 nodes unless stated otherwise. Each node 
corresponds to one neuron, whose dynamics is governed by 
the noise-driven Rulkov map. In Eq. (1) e l J = 1 if neu- 
ron i is coupled to neuron j and e i J = otherwise. Fol- 
lowing ref. 115611 . the preferential attachment is introduced via 
the probability II, which states that a new node will be con- 
nected to node i depending on its connectivity fc, according to 
n(fcj) = kjt as demonstrated schematically in Fig.Q] 

This growth and preferential attachment scheme yields a net- 

fc- 
work with an average degree k av = ^ ' , and a power-law 

degree distribution with the slope of the line equaling w — 3 on 




FIG. 1: Schematic presentation of growth and preferential attachment, as proposed by Barabasi and Albert 15611 . Each new node (white) 
preferentially attaches to two (thus here k av = 4) old nodes (black) that already have many other connections at that time. 



a double-logarithmic graph. Notably, analytical estimations 
predict the slope of the line to equal —3 II56I1 . We will use 
Barabasi-Albert scale-free networks having k av = 4 through- 
out this work (see Fig. [O, unless stated otherwise. 

In order to quantitatively study the degree of spatiotempo- 
ral synchronization in the network, and thus support below 
presented visual assessments of front propagation via space- 
time plots, we introduce, by means of the standard deviation, 
a synchronization parameter a (see e.g. [57]), which can be 
calculated effectively according to: 



1 

T ^ 
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a(n), CT (n) = -]>>>)]M-5>>)] 



i=l 



(2) 

It turns out that a is an excellent indicator for numerically 
measuring the spatiotemporal synchronization of excitations, 
hence revealing different synchronization levels and related 
transitions. From Eq. (2), it is evident that the smaller the 
synchronization parameter er, the more synchronous the neu- 
ronal network. Accordingly, when a — the network reaches 
complete synchrony. Final results shown below were aver- 
aged over 20 independent runs for each set of parameter val- 
ues to warrant appropriate statistical accuracy with respect to 
the scale-free network generation and numerical simulations. 



m. RESULTS 

We start by presenting space-time plots obtained with a 
fixed information transmission delay r = 700 and noise inten- 
sity w = 0.015, but different values of the coupling strength 
D. Results shown in the left three panels of Fig. [2] illustrate 
the spatiotemporal dynamics of neurons on the scale-free neu- 
ronal network having k av — 4. Evidently, for small coupling 
strengths [see panel (a), left] the excitatory fronts are quite 
nicely ordered in both time and space. However, both the tem- 



poral and spatial regularity increase further and substantially 
as D is increased [see panels (b) and (c), left]. Interestingly, 
by setting the information transmission delay to r = 1000 
and keeping the same noise intensity w = 0.015, the order in 
both time and space deteriorates substantially, as depicted by 
the right three panels of Fig. [2] Nevertheless, increasing the 
coupling strength can still improve the overall regularity of 
the excitatory fronts [comparing panels (a), (b) and (c), right]. 
It is thus revealed that different information transmission de- 
lays have a profound impact on the spatiotemporal regularity 
of excitatory fronts, whereas increasing the coupling strength 
always leads to an improvement of temporal and spatial syn- 
chronization. 

To investigate the impact of different information transmis- 
sion delays, outlined in Fig. [2] more precisely, we show in 
Fig. [3] space-time plots obtained with fixed coupling strength 
D = 0.01 and noise intensity w = 0.015, but different values 
of r. It can be observed that the spatiotemporal dynamics is 
ordered nicely if r — [see panel (a)]. When r = 200, this 
deteriorates drastically [see panel (b)], but is again revived at 
r = 600 [see panel (c)]. In fact, by closely examining space- 
time plots obtained with r = and r = 600, respectively, one 
can observe that the non-zero yet appropriately tuned informa- 
tion transmission delay can further enhance the regularity of 
excitatory fronts as compared to the case of r = 0. Quite re- 
markably, when r = 1000 the regularity of excitatory fronts is 
again heavily impaired [see panel (d)], yet with r — 1400 the 
order is restored anew [see panel (e)]. Indeed, the informa- 
tion transmission delay induced transitions to superbly syn- 
chronized neuronal activities on scale-free networks seem to 
appear intermittently, at roughly integer multiples of a given 
value of r, which equals approximately 600 — 700 in Fig. [3] In 
accordance with this preliminary assessment, it is expectable 
that with t = 1800 disorder in the temporal as well as the spa- 
tial domain sets in again [see panel (f)]. Visual investigations 
of Fig. 0thus reveal that regular and irregular front propaga- 
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FIG. 2: Left: Space-time plots of ar l '(n) obtained for r = 700 and w — 0.015 with different coupling strengths D, equaling: (a) 0.004, (b) 
0.008, (c) 0.016. Right: Space-time plots of x^\n) obtained for r = 1000 and w = 0.015 with different coupling strengths D, equaling: (a) 
0.004, (b) 0.008, (c) 0.016. In all panels, the system size is i = 1, 2 200 = TV. 
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FIG. 3: Space-time plots of (n) obtained for D = 0.01 and w = 
0.015 with different information transmission delays r, equalling: 
(a) 0, (b) 200, (c) 600, (d) 1000, (e) 1400, (f) 1800. In all panels the 
system size is i = 1,2..., 200 = N. 



tion appears intermittently as the delay is increased. Hence, 
it can be stated that finite (non-zero) information transmission 
delays play a pivotal role in the generation of spatiotemporal 
patterns of neuronal activity on scale-free networks. 

In what follows, the degree of spatiotemporal synchroniza- 
tion will be studied quantitatively via a (see Eq. 2), so as to 
support and validate the above visual assessments. Further- 
more, it remains of interest to examine the impact of different 
w and k av . In Fig.HJa), we first plot a in dependence on D for 
three different values of r. As visually interpreted by space- 
time plots presented in Fig. [2] larger coupling strengths indeed 
facilitate spatiotemporal synchronization in an monotonous 
manner. That is, as D increases, a decreases (irrespective 
of t), which is in agreement with previous studies examin- 
ing synchronization phenomena in neuronal as well as many 
other nonlinear systems. More eventful are results presented 
in Fig.|4jb), where a is presented in dependence on r for three 
different coupling strengths D. It can be observed clearly that 
certain values of r significantly facilitate spatiotemporal syn- 
chronization of excitatory fronts on scale-free neuronal net- 
works. The two minima of a appear at r w 700 and r w 1400, 
respectively, and are largely independent of D. This confirms 



the above claim that the information transmission delay in- 
duced transitions to spatiotemporally synchronized neuronal 
activity appear intermittently, at integer multiples of the given 
value of t. On the other hand, values of r outside these re- 
gions significantly impair synchronization, as can be inferred 
from the rather sharp ascends towards larger a beyond the op- 
timal delays. 

The delay-induced synchronization transitions, as well as 
the impact of increasing D, can be presented succinctly also 
via contour plots of u in dependence on the two parameters. 
Figure [5] features two such graphs obtained for w = 0.015 
[panel (a)] and w = 0.03 [panel (b)]. Transitions to highly 
synchronized states in dependence on the information tran- 
sition delay are clearly visible as extensive white regions 
(denoting smallest values of a) occurring at t k 700 and 
t k 1400, corresponding to the two minima depicted pre- 
viously in Fig. 5|b). Moreover, results presented in Fig. [5] 
clearly convey the impact of increasing the coupling strength 
D. In fact, not only do increasing D decrease a, as outlined 
above when interpreting Fig. Ufa), but also they broaden the 
span of r within which synchronous spatiotemporal neuronal 
activity is warranted. Notably, the optimal values of r shift 
insignificantly during the broadening. These features can be 
inferred also from Fig. HJb), yet the contour plots in Fig. [5] 
convey them more clearly. Finally, it is interesting to note that 
different values of w do not evoke significantly different re- 
sults, as can be appreciated by comparing panels (a) and (b) of 
Fig-HJ From this, we conclude that the delay-induced transi- 
tions to synchronous neuronal activity on scale-free networks 
are largely independent of D (apart from the broadening of the 
interval of r warranting optimal synchronization) and robust 
against reasonable variations of the noise intensity. 

Thus far, we have considered only scale-free networks hav- 
ing k av — 4. Since the average degree determines the density 
of interneuronal links, and is thus arguably an important pa- 
rameter, we present effects of different k av in Fig. [6] Panel 
(a) features space-time plots obtained with D = 0.004 and 
t = 500 for increasing average degree from left to right. Vi- 
sual inspection reveals that the impact of increasing k av is 
comparable to the impact of increasing D (comparing with 
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FIG. 4: (color online) (a) Dependence of the synchronization parameter a on D with different r. (b) Dependence of a on r with different D. 
Where applicable, other parameters are the same as in Fig. [2] 




FIG. 5: Contour plots of a in dependence on D and r with: (a) w = 0.015 and (b) w = 0.03. Information transmission delay induced 
synchronization transitions are clearly visible and largely robust to variations of the noise intensity w. 



the space-time plots presented in Fig. |2]i in that the excita- 
tory fronts propagate increasingly ordered in both time and 
space as the average degree increases. This can be confirmed 
quantitatively via a, as shown in Fig. |6jb). Indeed, larger 
k av shift lower the whole outlay of a in dependence on r, 
thus indicating improvement in the quality of spatiotemporal 
synchronization of neuronal activity on scale-free networks. 
Importantly, however, the oscillating outlay of a, along with 
the optimal value of r, is preserved and does not vary in 
dependence on k av . Therefore we conclude that the delay- 
induced synchronization transitions are a robust phenomenon 
of neuronal dynamics on scale-free networks, which an impor- 
tant role in achieving synchronized information transmission 
among neighboring as well as distant neurons. 

In order to further test the generality of our findings, we 
examine the impact of different N as well as different scal- 
ing exponents characterizing the scale-free degree distribu- 
tion. While the growth and preferential attachment algo- 
rithm proposed by Barabasi and Albert [56] yields a power- 
law degree distribution with the slope of the line w —3 on a 
double-logarithmic graph and serves as the most commonly 
used structure for testing theoretical models, also relevant for 



neuronal networks are slopes around —2, as reported in lfl8ll . 
Accordingly, we employ an alternative algorithm for scale- 
free network generation based on assigning a quenched fit- 
ness value to every node, and drawing links among them with 
a probability depending on the fitness of the two involved sites 
[15811 . Using exponentially distributed fitness and a threshold 
rule for linking nodes, we obtain a scale-free network with the 
scaling exponent —2, as presented in Fig. 3 of ifBUl . Results 
presented in Figs. [7f a ) an d (b) show clearly that variations of 
the system size do not notably influence the outcome of our 
simulations. In fact, the minima of a remain located at the 
same values of r and are of the same depth (with a reasonable 
error margin) irrespective of N. Likewise, by changing the 
scaling exponent from —3 (as given by the Barabasi-Albert 
algorithm) to —2 (as given by the algorithm described in 115 811 ) 
the results also do not change significantly in that the minima 
of a appear by roughly the same values of r [compare pan- 
els (a) and (b) of Fig. [7). Notable as a qualitative difference 
between the two scaling exponents characterizing the under- 
lying scale-free topology is the very first local minimum of 
<j by t ~ 250, occurring if the slope of the degree distribu- 
tion equals —2, as depicted in Fig. [7Jb). We conjecture that 
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FIG. 6: (color online) (a) Space-time plots of x^'(n) obtained with D = 0.004, w = 0.015 and r = 500 for different average degrees k a v, 
equalling (from left to right) 4, 8 and 12, respectively. In all presented panels the system size is i — 1,2 ... , 200 = N. (b) Dependence of a 
on t with different kav Other parameter values are the same as in panel (a). 




FIG. 7: (color online) (a) Dependence of a on r with different N of the Barabasi-Albert scale-free network l56ll having k av = 4. Other 
parameter values are D — 0.004 and w = 0.015. (b) Dependence of a on r with different N of the scale-free network generated as proposed 
in f58h (see also text for details). Other parameter values are D — 0.008 and w = 0.015. 



this minimum may be related to the subharmonic of the opti- 
mal t (1/3 of the first minimum at r « 740). Nevertheless, 
results presented in Fig. [7J attest to the fact that reported syn- 
chronization transitions on scale-free neuronal networks due 
to finite information transmission delays are largely robust to 
variations of the system size and of the scaling exponent char- 
acterizing the scale-free degree distribution. 

Finally in this section, we provide an explanation for the 
emergence of the newly reported delay-induced synchroniza- 
tion transitions. Up to now we have shown that the optimal 
value of r, resulting in the occurrence of the first minimum 
of a, as well as its reappearance at integer multiples, does not 
vary significantly in dependence on D [see Fig.HJb)], w (see 
Fig. [3]l, k av [see Fig. [6jb)], N [see Fig. [71 a)] or the scaling ex- 
ponent characterizing the underlying scale-free network [see 
Fig- 13b)]. This leads to the conclusion that parameters deter- 
mining the global neuronal dynamics do not play a significant 
role. Hence, one may adjust the local dynamics of each neu- 
ron by varying (3 and 7 (thus far we have not varied them). 
These two parameters affect the speed of the temporal evo- 
lution of yi(n), and consequently the predominant oscillation 
period of excitations. Results presented in Fig. [8] clearly show 



that the locations of minima of a shift to different values of r 
as (3 and 7 are varied. We thus determined the predominant 
oscillation period t osc of individual neurons within the scale- 
free network, and remarkably found that for j3 = 7 = 0.0006 
it is t osc w 1200, for (3 = 7 = 0.001 it is t osc rj 730, and for 
(3 = 7 = 0.0015 it is t osc ks 580. These data agree very well 
with the occurrence of the first minimum of a in dependence 
on t as presented in Fig. [8] Accordingly, we also conclude 
that the information transmission delay induced transitions to 
spatiotemporal synchronization of neuronal activity are due to 
the locking between r and the predominant oscillation period 
of individual neurons on the scale-free network. 



IV. SUMMARY AND DISCUSSION 

In sum, we have studied front propagation and synchroniza- 
tion transitions on scale-free neuronal networks in dependence 
on the information transmission delay, coupling strength and 
the average degree. We found that an appropriately adjusted 
delay length can significantly enhance synchronization of ex- 
citatory fronts in an intermittent fashion in dependence on r. 
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FIG. 8: (color online) Dependence of a on r for different com- 
binations of (5 and 7. Other parameter values are: D = 0.018, 
w = 0.015 and k av — 4 (Barabasi- Albert scale-free network 1 56]). 

The intermittent outlay emerges due to the locking between 
the delay and the inherent oscillation frequency of individual 
neurons on the scale-free network. Thus, approximately at ev- 
ery multiple of the inherent oscillation period of each neuron, 
the information transmission delay between coupled neurons 
results in supremely ordered and synchronized fronts of ex- 
citations. The widths of these dynamical regimes in depen- 
dence on r can be broadened, and the synchronization fur- 
ther improved, if the coupling strength or the average degree 
of the network is enlarged. In addition, we have examined 
the robustness of these findings to different levels of additive 
noise, as well as to different system sizes and scaling expo- 



nents characterizing the scale-free topology, finding that all 
the conclusions prevail. We have shown that fine-tuned infor- 
mation transmission delays can effectively supplement some 
recently identified mechanisms for the enhancement of syn- 
chronization O59ll60ll . as well as weak signals in general Il6lll . 
on scale-free networks. These conclusions seem to be sup- 
ported by real biological data stating that conduction veloci- 
ties along axons connecting neurons vary from 20 to 60 m/s 
dlH]. Real-life transmission delays are within the range of mil- 
liseconds, suggesting that substantially lower or higher values 
may be preclusive for optimal functioning of neuronal tissue. 
In future studies, it would be interesting to examine the impact 
of synaptic noise and different conductance states [63] on syn- 
chronization transitions in delayed complex networks, as well 
as to pinpoint the precise role of different aspects of structure 
and functioning of active neuronal networks 16411 . We hope 
that our present study will be a useful source of information 
when striving towards these goals. 
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